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ABSTRACT 

When  a  low  Rossby  number  barotropic  flow  accelerates  in  the  laterally  converging  half  of  a  strait,  the  local 
propagation  speed  of  long  topographic  waves  can  be  reduced  to  zero,  thereby  Mocking  or  preventing  the  formation 
of  a  steady  flow  downstream  from  the  strait.  An  inviscid  longwave  theory  is  presented  for  the  new  steady 
upstream  and  downstream  states  that  evolve  from  the  blocking  wave.  The  enhanced  inshore  cyclonic  vorticity 
extending  far  downstream  suggests  that  topographic  jetogenesis,  rather  than  lateral  eddy  diffusion,  in  major 
ocean  straits  (e.g.,  Yucatan  and  Florida)  may  be  important  in  generating  or  reforming  boundary  currents. 


I.  Introduction  and  simple  considerations 

Steady  barotropic  currents  tend  to  follow  continental 
isobaths  when  the  curvature  and  the  Rossby  number 
R  are  sufficiently  small.  In  a  region  of  converging  iso¬ 
baths  mass  conservation  requires  the  downstream  ve¬ 
locity  to  increase,  and  potential  vorticity  conservation 
requires  cross-isobath  motions  to  account  for  the  rel¬ 
ative  vorticity  changes.  Thus  the  velocity  field  in  a  con¬ 
verging  channel  can  increase  to  a  critical  state  such 
that  the  local  propagation  speed  of  long  topographic 
waves  vanishes,  in  which  case  no  steady  downstream 
state  exists  for  the  given  upstream  flow.  A  theory  will 
be  developed  for  the  new  upstream  and  downstream 
state  as  a  function  of  the  topographic  convergence. 

The  “blocking”  effect  that  occurs  when  the  flow  be¬ 
comes  topographically  critical  is  analogous  to  that  in 
classical,  nonrotating  open  channel  hydraulics  ( Long 
1972)  when  an  obstacle  is  towed  along  the  bottom  of 
an  otherwise  resting  layer  of  water  with  a  speed  greater 
than  that  of  long  surface  gravity  waves.  A  wave  of  ele¬ 
vation  then  forms  with  its  trailing  edge  fixed  near  the 
obstacle,  and  with  its  leading  edge  propagating  ahead 
of  the  obstacle.  In  this  way  a  new  upstream  state  with 
an  increased  water  level  is  produced,  as  well  as  a  new 
( supercritical )  downstream  state.  Similar  effects  occur 
in  a  rotating  stratified  fluid  when  it  is  blocked,  and  here 
it  is  the  Kelvin  wave  propagation  that  is  relevant.  {See 
the  reduced  gravity  model  of  Gill  ( 1977 ).] 

The  control  effect  of  continental  topography  on  a 
parallel  current  has  been  studied  by  several  authors. 
Hughes  ( 1985,  1986,  1987)  has  computed  particular 
“conjugate”  flow  states,  for  stratified  as  well  as  baro- 
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tropic  flow,  having  the  same  specified  potential  vorticity 
on  the  same  streamline.  Hydraulic  transitions  between 
such  end  states  are  possible,  but  the  theory  does  not 
address  the  question  of  the  new  upstream  and  down¬ 
stream  states  that  evolve  when  an  arbitrary  upstream 
flow  is  blocked.  In  that  case  the  spatial  distribution  of 
an  entire  field  of  potential  vorticity  is  modified,  so  that 
the  problem  becomes  more  difficult  than  the  classical 
one  (Long  1972)  in  which  only  one  or  two  parameters 
(e.g.,  the  fluid  height)  are  necessary  to  describe  the 
upstream  state.  For  this  reason  it  is  desirable  to  restrict 
our  problem  to  very  simple  fields  of  upstream  potential 
vorticity,  such  that  only  a  finite  number  of  degrees  of 
freedom  are  involved  in  its  modification.  Herein  lies 
the  importance  of  p’ecewise  uniform  potential  vorticity 
flow  in  the  following  inviscid,  barotropic.  and  /-plane 
model. 

Let  us  start  by  considering  the  steady  flow  shown  in 
Fig.  1 ,  in  which  the  distinctive  feature  of  the  topography 
is  the  relatively  large  convergence  of  the  isobaths  near 
the  central  axis  of  the  channel,  as  compared  with  the 
convergence  near  the  sidewalls.  The  downstream  (x* ) 
gradient  of  isobathic  depth  h*(  x*.  y* )  is  assumed  smali 
compared  to  the  cross-isobath  (>•*)  gradient,  so  that 
—du*j  dy*  may  be  used  as  the  long-wave  approxi¬ 
mation  for  the  relative  vorticity,  where  «*  is  the  di¬ 
mensional  downstream  velocity,  and  v*  <s  u*  is  the 
y*  velocity  component.  Far  upstream  in  Fig.  I  a  uni¬ 
form  u*  =  U*  is  assumed,  with  Rossby  number 
UtlfW*  <3  1 ,  where  />  0  is  the  Coriolis  parameter. 
Under  this  condition  a  first  approximation  to  the  law 
of  conservation  of  potential  vorticity  requires  the  steady 
streamlines  to  coincide  with  the  isobaths,  as  previously 
mentioned.  To  conserve  mass  between  laterally  con¬ 
verging  isobaths  the  velocity  U*  on  the  central  axis 
must  increase  relative  to  U* ,  whereas  the  wall  velocities 
U* ,  U*  change  less.  This  implies  that  a  jet  forms  in 
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the  converging  half  of  the  strait,  with  the  necessary 
relative  vorticities  being  produced  by  cross-isobath 
stretching  or  squashing  of  columns.  Thus  the  cyclonic 
vorticity  at  P,  (Fig.  1 )  and  the  anticyclonic  vorticity 
at  P2,  require  cross-slope  displacements  as  indicated 
by  the  arrows. 

All  but  one  of  the  following  calculations  will  be  done 
for  models  (Fig.  2)  in  which  the  smooth  continental 
slopes  are  replaced  by  escarpments  (the  dashed  curves), 
across  which  the  nondimensional  depth  changes  dis- 
continuously  from  h  =  1  to  h  -  &  >  1.  The  lateral 
length  scale  is  some  W*,  chosen  (as  convenient)  to 
make  either  W0  or  W  j  or  W2  equal  to  unity  in  Fig.  2, 
and  the  nondimensional  D(x)  gives  the  downstream 
(x)  variation  in  width  of  the  deepest  part  of  the  channel 
(the  distance  between  escarpments).  The  uniform 
upstream  speed  U *  is  the  velocity  scale,  L’t/W*  is 
the  vorticity  scale,  and  the  Rossby  number  is  R 
=  Ut/fW*. 

The  main  problem  is  to  predict  the  steady  flow  hav¬ 
ing  three  piecewise  uniform  potential  vorticity  layers, 
and  a  given  total  volume  flux  at  each  jc.  These  three 
layers  are  separated  by  two  interfaces  whose  displace¬ 
ments  L i,  L7  from  their  respective  escarpments  are 
tentatively  assumed  positive  in  the  —y  direction  (to 
anticipate  the  jetlike  profile  in  Fig.  2). 

Since  the  upstream  relative  vorticity  vanishes,  and 
since  potential  vorticity  is  conserved,  the  cross-escarp¬ 
ment  stretching  (£1  -  1)  of  columns  in  the  L,  layer 
produce  the  nondimensional  vorticity 


Fig  .  1 .  Schematic  diagram  for  steady,  subcritical  ( low  upstream 
Rossby  number)  barotropic  flow  through  a  strait.  The  sketched  iso¬ 
baths  converge  more  rapidly  in  midchannel  than  the  isobaths  near 
the  shores.  The  “zero-order"  tendency  for  the  flow  to  follow  the 
isobaths  suggests  a  large  U*  must  form  in  the  highly  convergent 
central  portion  of  the  straits.  The  implied  relative  vorticities  require 
cross-isobath  displacements  at  points  P,  and  P2  as  shown  by  the  ar¬ 
rows;  W*  is  a  characteristic  channel  width  (dimensional). 


u*»i 


Fig.  2.  An  escarpment  model  as  a  limit  of  Fig.  I.  The  depth  is 
uniform  on  the  two  "shelves"  of  uniform  nondimensional  widths 
H'|,  H'2.  The  nondimensional  depth  increases  (discontinuously)  from 
unity  to  d  across  the  escarpments  (dashed).  Cyclonic  vorticity  f  is 
generated  in  the  strip  of  width  L,  located  between  the  upper  escarp¬ 
ment  and  the  interface  (solid  curve)  separating  the  layers  of  uniform 
potential  vorticity.  Anticyclonic  vorticity  is  generated  in  the  L2  layer. 


Likewise,  the  magnitude  of  the  anticyclonic  vorticity 
in  the  L2  layer  is  f_  =  ( 1  -  (31 )//?,  or 

(1.2) 

According  to  the  long-wave  approximation  the  down¬ 
stream  velocity  varies  linearly  in  these  two  vortical  lay¬ 
ers,  and  outside  the  velocities  Lj,  £/,  and  U2  are  in¬ 
dependent  of  y,  where  V  is  the  uniform  velocity  in  the 
central  region,  U\  is  the  upper-shelf  velocity,  and  U2  is 
the  lower-shelf  velocity. 

It  is  instructive  to  consider  first  the  simplest  example, 
in  which  both  IF, ,  W2  are  much  larger  than  =  1 , 
in  which  case  mass  conservation  in  the  relatively  wide 
regions  (W,,  W2)  requires  V\  =  U2  I.  Adding  this 
to  the  appropriate  shear  [Eq.  ( 1 . 1 )  or  ( 1 .2)]  gives  the 
velocity  U  between  the  two  interfaces 

1  +  fL,  =  U=  1  +L2f_,  (1.3) 

and  then  Eq.  (1.2)  yields 

J-2  ~  0L] .  (1.4) - 

At  any  section  x  the  total  volume  transport  between 
the  two  interfaces  is  the  sum  of  two  parts,  one  of  which 
[0 U(D  -  L| )]  occurs  in  deep  water,  and  the  other 
[  |  ( U  +  1  )L2\  occurs  in  the  shallow  anticyclonic  layer. 

TTiis  sum  is  independent  of  jc  and  equal  to  the  transport - 

/3W0  =  0  at  x  =  -oo,  so  that  ~ 

1  =(1  +  {£,)(/>- L,)  +  (l  +[L1/2)L2/P.  ~ 

The  use  of  ( 1 .4 )  then  gives  the  quadratic  equation  ’  odes 
L]  —  2DLi  +  2(1-  D)/i  =  0  (1.5)  or 
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and  the  solution  satisfying  =  0  when  D  -  1  is 

L,  =  T>-(D:--2<1  -  D)/f),/J.  (1.6) 


This  is  positive  because  D(x)  <  1.  and  thus  our  initial 
assumption  Lt  >  0.  L2  >  0  is  verified.  The  discriminant 
of  the  quadratic  vanishes  at  a  section  where 


D(x) 


r  +  U2 


2 

-f — 


1/2 


<  1, 


in  which  case  L,  =  D,  L2  =  t$D,  and  consequently  a 
triangular  jet  forms  having  a  maximum  velocity 


U=  1  +  fL,  =  (1  +  2ft1'2  = 


1/2 


located  on  the  lower  escarpment.  If  R  =  U*  / /  W  * 

1 ,  where  the  central  channel  width  W  *  now  provides 
the  tV*  scale,  this  U  is  very  large  compared  to  the 
upstream  current,  and  thus  we  see  that  pronounced 
in  viscid  jetogenesis  can  occur  in  the  laterally  converging 
half  of  a  channel.  The  large  values  of  the  cyclonic  and 
anticyclonic  vorticities  may  then  cause  separation  as 
the  jet  enters  the  diverging  half  of  the  channel  (Stern 
and  Whitehead  1990).  This  is  one  kind  of  transition 
that  would  prevent  a  reversal  of  the  vorticity  generation 
in  the  diverging  half  of  the  strait,  and  that  would  allow 
a  jet  to  appear  far  downstream. 

For  the  smaller  tV, ,  W2  ( narrower  shelves  in  sections 
2  and  3),  a  hydraulic  transition  or  blocking  can  occur 
first,  and  the  single  shelf  (W2  =  0)  model  in  section  2 
will  address  the  main  question,  namely,  what  happens 
if  min D(x)  is  small  enough  to  block  the  uniform  up¬ 
stream  current. 

Both  kinds  of  transition  seem  to  be  relevant  to  the 
influence  of  major  ocean  straits  on  the  circulation  in 
the  basins  they  connect.  One  example  is  the  Atlantic 
western  boundary  current,  which  appears  to  become 
very  diffuse  (and  remote  from  boundaries)  after  en¬ 
tering  the  Caribbean  Sea,  and  until  it  reforms  as  a  jet 
in  the  vicinity  of  the  Straits  of  Yucatan  (Gordon  1967; 
Molinari  et  al.  1981 ).  This  Gulf  Stream  jet  then  passes 
through  the  Gulf  of  Mexico,  and  enters  the  Straits  of 
Florida  where  the  maximum  mean  inshore  vorticity 
increases  to  a  value  approximately  equal  to  the  Coriolis 
parameter  ( Brooks  and  Niiler  1977 ).  It  is  unlikely  that 
this  large  inshore  vorticity  can  be  correctly  explained 
by  a  lateral  diffusion  mechanism  and  a  proper  expla¬ 
nation  should  be  based  on  an  inertial  (potential  vor¬ 
ticity  conserving)  mechanism,  taking  into  account  the 
isobaths  (and  the  isopycnals)  in  the  strait  (Stern  1991 ). 


2.  A  converging  channel  with  a  single  escarpment 

In  Fig.  3a  tV2  has  been  set  equal  to  zero,  W0  =  1  is 
the  length  unit,  and  D(x)  1  is  the  distance  between 
the  curved  lower  wall  and  the  escarpment.  Note  that 
in  a  long-wave  theory  the  curvature  of  isobaths  '«  un¬ 
important,  compared  to  their  y  separation,  so  that  Fig. 


3a  is  dynamically  equhaient  to  Fig.  2  when  tt'2  -  0. 
(Such  a  model  might  be  realized  in  a  long  rotating 
channel  by  towing  an  obstacle  placed  on  a  wall  through 
an  otherwise  resting  liquid.)  In  this  case  a  single  inter¬ 
face  separates  two  regions  of  uniform  potential  vortic¬ 
ity,  and  R  =  V*  /  j  H  *  and  fi  are  first  assumed  to  be 
such  that  a  steady  long-wave  solution  exists  with  /. 
-  L(D(x))  >  0.  As  previously  indicated  ( 1. 1 ),  the  rel¬ 
ative  vorticity  immediately  above  this  interface  is 

f  >  0  (2.1) 

A 

and  below  the  interface  the  velocity  V  is  independent 
of  y.  The  local  volume  transport  U(D  -  L)0  in  the 
latter  region  must  equal  the  value  (0)alx  -  -  % ,  and 
therefore 


V  = 


1 

D  -  L  ' 


(2.2) 


The  shear  above  the  interface  ( Fig.  3a)  gives  the  velocity 

Ut  =  U-a  (2.3) 

at  and  above  the  escarpment.  Therefore  the  total  vol¬ 
ume  transport  above  this  interface  is  U,  W  j  +  1  ( Lj 
+  U)Lfi,  and  the  latter  must  equal  the  transport  H  j 
atx=  -oo.  Simplification  of  this  equality  using(2.2)~ 
(2.3)  yields 


1  = 


1  +  LpJWj 
D- L 


-  r l 


2  Wi 


(2.4a) 


D  -  L  + 


1  +  Iff  l  FFj 
1  +  XL  +  /?T2f/2Wj  ' 


(2.4b) 


The  implicit  solution  (2.4b)  for  L  =  L(D)  reduces  to 
D  =  1  +  L((5jWx  —  £+!)+•■•  when  L  -*•  0,  and 
since  D  <  la  necessary  condition  for  the  L  >  0  as¬ 
sumption  is 


or 


f>  I  +0/W',  (2.5a) 


R  < 


0~  1 

1  +  /3/H  j  ■ 


(2.5b) 


The  right-hand  side  of  ( 2.4b)  is  positive  and  asymptotes 
to  D  =  L  when  L  is  large,  so  that  a  minimum  D  =  D, 
exists.  This  point  separates  two  branches  ( Fig.  3b)  hav¬ 
ing  different  L  values  for  the  same  D,  and  if  mini) 
<  Dc  no  steady  solution  exists  for  the  upstream  flow 
with  uniform  velocity. 

For  relatively  large  H’1  Eq.  u.4aj  reduces  to  "  qua¬ 
dratic  equation 


f 2  - 


t(n-L\ 

\  a 


1  i  i  --  o)n  -  o 


(2.6) 
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<  a  ) 


The  minimum  width  Dc  of  the  strait  necessary  for 
blocking  in  this  case  is  given  by  setting  the  discriminant 
equal  to  zero,  that  is. 


(c-e)! 


4(1  -  A)/f  -  I  A  + 


1\2  4 


n  r 


and  therefore 


For  large  W ,  Eqs.  (2.2)-(2.4a)  also  imply 


Fig.  3.  (a)  Flow  ofa  uniform  upstream  current  U0  =  1  in  a  channel 
with  a  single  escarpment  (dashed  horizontal  line).  The  distance  D(x) 
of  this  from  the  lower  vertical  wall  varies  slowly  on  the  W0  =  1  length 
scale.  The  interface  between  the  piecewise-uniform  potential  vorticity 
layers  is  displaced  downward  (negative  y)  by  an  amount  L(x)  >  0 
in  the  steady  state.  Above  the  interface  nondimensional  vorticity  f 
is  thereby  generated,  and  the  velocity  varies  linearly  between  the 
values  U\  and  U.  (b)  Sketch  of  a  long  wave  solution  for  L  as  a 
function  of  D.  The  branch  point  at  D  =  Dc  is  the  smallest  value 
possible  foragiven  Rossby  number  at  x  =  -oo  .(c)  Schematic  diagram 
of  the  temporal  evolution  of  the  potential  vorticity  interface  ( L(x , 
0)  =  0]  when  a  small  amplitude  constriction  ( the  stipled  “strait”)  is 
suddenly  moved  leftward  at  t  =  0  with  a  uniform  subcritical  speed. 
The  transverse  velocities  v  forced  at  /  =  0  tend  to  deflect  the  interface 
to  L(  x,  1 1 ),  and  this  generates  relative  vorticities  whose  sense  is  in¬ 
dicated  by  the  circular  arrows.  These  vortices  cause  the  interface  to 
propagate  to  L(x ,  t2)  for  /2  >  /, .  As  t  -»  *  the  ridge  ahead  of  the 
constriction  continues  propagating  to  x  =  -oo,  whereas  the  trough 
remains  behind  in  equilibrium  with  the  v  forced  by  the  constriction, 
as  shown  in  (a). 


and  the  solution  that  vanishes  at  D  =  1  is 


£/,-*•  1,  (2.9a) 

which  is  obviously  required  to  satisfy  the  transport  re¬ 
lation  W{U\  =  W\  in  the  wide  1 )  region. 

The  question  arises  as  to  the  physical  consistency  of 
an  inviscid  escarpment  model,  and  its  relation  to  a 
model  with  slightly  smoother  topography.  Accordingly, 
the  critical  value  for  (2.4b)  has  been  compared  (in  the 
Appendix)  with  the  corresponding  value  for  a  model 
having  a  smoother  continental  slope.  The  numerical 
results  indicate  that  the  inviscid  escarpment  limit  is 
approached  asymptotically  for  large  finite  continental 
slopes  (in  which  case  the  lateral  viscous  stresses  might 
be  consistently  neglected.)  This  justification  of  the  es¬ 
carpment  model  is  necessary  because  it  provides  the 
only  tractable  nonlinear  theory  for  obtaining  the  new 
upstream  state  when  the  old  one  is  blocked,  as  will  be 
seen.  Before  proceeding  to  this  question,  a  physical  ex¬ 
planation  of  Fig.  3a,  and  motivating  arguments  for  Fig. 
4  are  in  order. 

The  appearance  of  a  steady-state  “trough”  (i.e.,  L 
>  0  in  Fig.  3a)  can  be  qualitatively  explained  in  terms 
of  well-known  (e.g.,  Stem  1 99 1 )  properties  of  free  and 
forced  topographic  long  waves.  Consider  the  evolution 
from  time  t  -  0  when  the  fluid  is  resting  with  L(x,  0) 
=  0,  and  then  the  strait  (constriction )  is  moved  leftward 
with  uniform  speed  (Fig.  3c ).  This  obviously  forces  an 
initial  transverse  velocity  v  >  0  at  x  <  0,  and  v  <  0  at 
x  >  0,  so  that  the  potential  vorticity  interface  at  small 
U  >  0  is  deflected  as  shown.  The  vorticity  generated 
by  the  cross-escarpment  displacement  then  causes  rel¬ 
ative  leftward  phase  (and  group)  propagation  of  L  if 
the  relative  upstream  velocity  is  subcritical  (i.e.,  smaller 
than  the  propagation  speed  of  free  topographic  long 
waves).  Therefore,  at  1 2  >  t\  the  L-ridge  moves  ahead 
of  the  strait  and  continues  propagating  to  x  =  — co  as 
an  unforced  (free)  wave.  But  the  advance  of  the  trough 
eventually  ceases  under  the  influence  of  the  forced  v, 
and  thus  the  steady-state  trough  L  >  0  (Fig.  3a)  is 
attained  in  the  strait. 

This  picture  applies  to  small  values  of  ( 1  -  Dm), 
and  for  larger  ones  a  first-order  nonlinear  hyperbolic 
equation  (Stern  !°91)  applies  to  the  evolutionary 
problem.  In  this  case  the  current  in  the  converging  strait 
may  increase  to  a  value  such  that  the  “local  L-propa- 
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gation  speed”  (according  to  the  method  of  character¬ 
istics)  vanishes,  at  which  point  (x)  the  trailing  edge  of 
the  ridge  remains  behind  in  the  vicinity  of  the  strait, 
while  the  leading  edge  continues  propagating  to  x 
=  -  qo  in  the  locally  subcritical  region  (D  1 ) .  Figure 
4  is  a  schematic  diagram  of  this  conjectured  flow,  as 
seen  by  an  observer  on  the  strait.  As  /  -*■  oo  the  up¬ 
stream  propagating  blocking  ridge  (A  >  0)  produces  a 
new  mean  field  with  an  anticyclonic  shear  layer  at  all 
finite  x  <  0.  Farther  downstream  a  transition  to  the 
supercritical  branch  (cf.  Fig.  3b)  of  the  new  steady  so¬ 
lution  occurs,  with  cyclonic  vorticity  appearing  in  the 
trough  region  (Fig.  4).  This  (incomplete)  picture  of 
the  temporal  evolution  provides  the  motivation  for  the 
following  “self-consistent”  steady  long-wave  theory. 

Let  Xo  in  Fig.  4  denote  the  ( unknown)  zero  crossing 
point  of  the  interface,  let  D0  denote  its  distance  from 
the  curved  wall,  and  let  Ao  denote  the  ( unknown )  hor¬ 
izontal  velocity  at  jco-  This  velocity  must  have  the  same 
value  in  both  layers  because  L(D0 )  =  0  implies  that 
there  is  no  vorticity  at  xq.  Since  the  total  volume  trans¬ 
port  between  the  two  walls  must  have  the  same  value 
at  all  t,  x  (even  at  x  =  ±co  in  the  transient  region ),  it 
follows  that  the  transport  To  Aid  +  U'0WX  at  x0  must 
equal  the  original  transport  Wx  +  /3  at  x  =  -oo ,  or 


occur  where  dD/dx  =  0,  that  is.  at  a  =  0  where  D  ~  D„, 
(Fig.  4),  or 


This  equation  determines  the  fin-4  unknown  in  the  1. 
=  L(  D)  solution,  and  gives  the  new  downstream  ve¬ 
locity  field.  With  D0  and  to  known,  the  mass  balance 
for  x  <  A'o  gives  A (D)  and  the  new  vorticity  field  up¬ 
stream  from  the  strait. 

This  outline  of  the  calculation  will  now  be  carried 
out  for  the  simplest  case  of  large  W',  (cf.  Eq.  (2.7 )].  in 
which  case  the  upper-layer  velocity  (2.9a)  is  nearly 
equal  to  unity  ( A0  =  1 )  at  all  x,  and  A  =  1  +  {"L 
beneath  the  interface  in  the  cyclonic  (x  >  x<>)  region. 
The  steady  mass  balance  for  this  region  then  becomes 
(1  +  ft)(Z)  -  L)  =  A>,  or 

L2-t|D-ij  +  (Ao-A)/f=0.  (2.12) 

For  the  (/?,/?)  values  listed  subsequently,  the  zero¬ 
crossing  point  is  at  A®  <  0  and  therefore  D  <  Do  in  Ay 
<  x  <  0.  Since  L  >  0  as  D  -  D0  -*■  0",  Eq.  (2. 12) 
requires 


l+P/Wi 

1  +  0Do/lVt  1 


(2.10) 


Let  L(x)  2*  0  be  the  interface  displacement  in  the  x 
>  Xq  region,  and  let  A(jc)  >  0  be  the  interface  displace¬ 
ment  in  the  x  <  Xo  region.  By  applying  mass  conser¬ 
vation  to  the  former  region  (as  done  previously),  and 
using  L(  A0)  =  0,  a  solution  for  L(D)  is  obtained  con¬ 
taining  one  unknown  (Ao).  Since  the  blocking  con¬ 
dition  min/)  <  Dc  is  satisfied,  the  solution  requires  a 
branch  point  with  dL/dD  =  co,  and  with  dL/dx  =  (dL/ 
dD)(dD/dx)  finite.  Therefore  the  branch  point  must 


Fig  .  4.  The  modification  of  the  upstream  flow  by  a  topographic 
blocking  wave.  The  minimum  width  Dm  of  the  deep  channel  is  at  x 
-  0.  A  wiggly  arrow  indicates  the  upstream  propagation  of  the  an¬ 
ticyclonic  vorticity  ( >  0)  layer  of  thickness  A(x)  >  0.  This  produces 
a  reduction  of  the  lower-layer  velocity  ( l/_)  relative  to  the  undisturbed 
velocity  { =  I )  at  x  =  -  o& .  At  jc  >  a®  a  compensating  cyclonic  vorticity 
(f  >  0)  of  thickness  L(x)  >  0  develops  downstream.  At  large  finite 
distances  upstream  from  the  obstacle,  and  at  large  finite  time  after 
startup,  a  steady  state  is  assumed,  even  tnough  the  blocking  wave 
continues  propagating  to  x  -*•  ~oo . 


1  R 

0>f  “FT- 


(2.13) 


and  therefore,  the  solution  that  vanishes  at  Xq  is 


=K) 

-M- 


4(A) -A)/f  >0.  (2.14) 


Downstream  from  the  branch  point  (x  =  0)  defined 
by  ( 2. 1 1 )  the  appropriate  solution  of  ( 2. 1 2 )  is 

2L  =  |/)-^+  \d-  j^-4(A)-Z))/f  ' 


and  the  branch  point  (D  =  Dm)  occurs  where  the  dis¬ 
criminant  vanishes,  or 

0  =  (Am-|J2-4(A)-  Dm)n 


Therefore  the  zero-crossing  point  of  the  interface  occurs 
where 


A,  =  ^(/)m  +  |)\  (2.16) 

Equations  (2.14)— (2.16)  give  the  new  downstream 
state  (L  >  0)  after  blocking  occurs;  but  (2.13)  restricts 
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the  solution  to  the  parametric  regime  1  <  ff>()  =  ( £Dm/ 
2  +  1/2)*,  or  iDm  >  1.  and  by  combining  this  with 
( 2.9 )  we  get 

r'  <Dm<2r,/2~i  '  (2-17) 


for  the  restricted  Dm  range. 

The  A  >0  solution  for  x  *£  .v«  is  obtained  by  noting 
that  V  =  1  -  Af.  is  the  velocity  below  the  escarpment 
in  big.  4,  where  -  f//?.  Therefore  the  total  transport 
below  the  interface  is 


If  D,„  is  less  than  the  lower  bound  in  ( 2. 1 7 )  then  a  sign 
change  in  .v0  is  suggested  (by  the  above  value  of  x0 
=  0 ),  in  w-hich  case  a  separate  calculation  is  necessary. 
The  foregoing  calculation  also  does  not  apply  when  R 
exceeds  the  bound  in  (2.5b),  because  the  current  at  x 
=  -  x  is  strong  enough  to  prevent  upstream  topo¬ 
graphic  wave  propagation. 

This  steady  longwave  theory  neglects  the  possibility 
of  persistent  short  wave  effects,  such  as  are  bound  to 
occur  (as  “shocks”)  in  the  temporal  evolution.  These 
eddy  effects  might  be  studied  in  a  finite  difference  com¬ 
putation  using  a  smooth  continental  slope  model,  such 
as  appears  in  the  Appendix. 


DpU-  +  -  ( 1  +  C/_)A  =  D03 


or 


20 


(D  -  Do)  —  0,  (2.18) 


1  >  D>  D0>  Dm>~.  (2.19) 

These  inequalities  imply  fl1  -  0D  <  0,  and  therefore 
the  solution  of  (2.18)  satisfying  A  0  is 


A  =  -  -  0D  + 


KM 


+  2 0(D  -  A))/f- 


-11/2 


(2.20) 

At  x  <  x0,  U-  -  1  -  Af_  becomes 

V-  =  D£  -  (f2D2  -  2CA>  +  1 )l/2,  (2.21) 

and  far  upstream  ( D  -*  1 )  the  velocity  is 

f/-(-oo)  =  C~(C2-2CO0+  1)1/2,  (2.22) 


where  D0  is  given  by  (2.16). 

These  results  may  be  summarized  as  follows.  No 
blocking  occurs  when  the  minimum  channel  width  Dm 
exceeds  the  upper  bound  in  (2.17),  and  in  this  case 
the  original  upstream  flow  merely  produces  a  sym¬ 
metrical  trough  L(x)>  0  in  the  vicinity  of  the  obstacle. 
When  Dm  is  slightly  less  than  the  upper  bound  in 
(2.17),  Eq.  (2.16)  gives  Do  T  which  means  x0  -*■ 
-x,  A  0  +  ,  and  L(x)  >  0  at  all  x.  But  there  is  a 
branch  point  at  x  -  0,  and  a  transition  from  (2.14)  to 
(2.15),  so  that  at  x  =  +<x  we  get  2  jL  =  2(1  -  I/O 
>  0.  If  Dm  is  decreased  further  below  the  upper  bound 
in  (2.17)  then  Z>0  <  1  in  (2.16),  x0  becomes  finite 
negative,  A(x)  >  0,  and  a  finite  modification  of  the 
upstream  state  (Fig.  4)  occurs  in  which  anticyclonic 
vorticity  appears,  in  addition  to  the  cyclonic  vorticity 
in  the  L  >  0  strip  which  extends  far  downstream.  When 
Dm  equals  the  lower  bound  (O')  in  (2.17)  then  Do 
-  1  /  C  =  Dm  so  that  x0  ~  0,  and  the  modified  upstream 
velocity  (2.22)  is 


u.  =  C-(C2  -  i),/2  <  1. 


3.  The  two-escarpment  model  of  a  strait 

Figure  2  has  a  greater  geometric  resemblance  to  a 
strait,  and  a  jetlike  (L,  >  0,  L2  >  0)  solution  is  sought 
when  W,  -  1  =  W2  provide  the  length  scale.  The  cy¬ 
clonic  vorticity  f  is  given  by  ( 1 . 1 ) ,  and  the  ( magnitude 
of  the)  anticyclonic  vorticity  is  given  by  (1.2).  The 
uniform  velocities  on  the  upper  and  lower  shelf  are 
given,  respectively,  by  U\  -  U  -  fZ-i,  U2  =  U  - 
where  U  is  the  uniform  velocity  between  the  interfaces. 
What  is  min(ZA)  for  a  hydraulic  or  for  a  separating 
transition  in  this  model? 

For  the  lowest  layer  (Fig.  2)  of  uniform  potential 
vorticity  and  unit  depth  the  total  volume  transport  ( U 
-  L2$- )( l  -  L2]  at  any  D(.x)  must  equal  unity,  or 

V  =  - — +  L2H- .  (3.1) 

1  -  L2 

In  the  middle  layer  between  the  two  interfaces,  the  total 
volume  transport  0U{D  -  Lt)  +  (U  —  L2H-/2)L2 
equals  0Wn,  or 

L^V  -  UD  -  Wo  +  /r'(l/~  L2$-/2)L2.  (3.2) 
In  the  uppermost  layer  the  transport  balance  is 

V-  U,  +0UJ-  r/-./2)L,  =  l.  (3.3) 

The  values  of  (/?,  R,  H  o)  necessary  to  justify  the  Lt 
>  0,  L2  >  0  assumpt’on  are  obtained  by  letting  x  ~* 
-00,  L,  —  0,  L2  -*  0,  U  ~  1  «  y  ~+  0,  B0  -  D 
=  2 M(x)  -*■  0  +  .  In  this  region  the  linearization  of 
(3.1 )— ( 3.3 )  yields 

g  =  (l  +  f-)T2, 

L,  =  nW0  -  2M  +  0~'t2. 

M  =  (f-J8)L„ 
and  consequently 

_ 2  _ 

2  1  +  r-  -tr-0)  id  +  f_)iFo  +  rl] ' 

(3.4) 
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The  last  equation  requires 


and  since  M  >  0  Eq.  (3.4)  requires 

<r-0)l(»  +  +  +  f-)>0.  (3.6) 

When  this  is  multiplied  by  /3  and  when  =  f  is  used, 
we  get 


(f-0)[(f+0)W'o  +  l]-(f  +  /8)>0 


or 


+  i£<f2-(j8_  0 

WQ  f 


R‘ 


(3.7) 


This  condition  for  positive  (L\ ,  L2),  as  well  as  the  one 
in  (3.5)  is  obviously  satisfied  for  sufficiently  small 
Rossby  number. 

The  significance  of  ( 3.7 )  appears  by  considering  the 
conditions  under  which  the  basic  flow  at  jc  =  —  oc  can 
support  a  free  (i.e.,  M(x)  =  0)  stationary  long  wave 
of  infinitesimal  amplitude.  Equation  ( 3.4)  implies  that 
this  will  occur  when  the  denominator  vanishes,  or  when 
the  left-hand  side  of  ( 3.6)  vanishes,  or  when  an  equality 
sign  replaces  the  inequality  in  (3.7).  The  latter  equality 
then  gives  the  Rossby  number  for  a  stationary  long 
wave,  and  for  smaller  R  the  long  wave  propagates  up¬ 
stream.  It  is  under  the  latter  condition  (3.7)  that  a 
steady  forced  solution  with  L,  >  0,  L2  >  0  for  D(x) 
<  W0  is  obtained.  (See  the  previous  argument  in  con¬ 
nection  with  Fig.  3c.) 


Table  la.  Critical  values  of  the  parameters  for  the  two  escarpment 
model  (section  3  and  Fig.  2)  when  ft  =  2.  and  for  various  values  of 
K,  U'o,  The  first  row  for  each  entry  is  the  minimum  D,  and  the 
following  rows  are  U,,  V,  U2,  respectively.  The  symbol  H  indicates 
a  hydraulic  or  blocking  transition. 


R 

Wo 

1 

2 

4 

8 

32 

1/4 

.85  H 

1.40  H 

2.5  H 

4.68  H 

26.08  H 

.544 

.168 

,178 

.170 

.769 

1.455 

1.831 

1.821 

1.829 

1.230 

1.167 

1.332 

1.327 

1.331 

1.080 

1/6 

.68  H 

1.16  H 

2.12  H 

4.36  H 

15.68  H 

.383 

.302 

.294 

.426 

.303 

1.961 

2.068 

2.078 

1.903 

2.066 

1.288 

1.327 

1.331 

1.268 

1.326 

1/10 

.52  H 

.92  H 

1.72  H 

3.30  H 

12.8  H 

.386 

.389 

.389 

.380 

.377 

2.507 

2.501 

2.500 

2.517 

2.524 

1.313 

1.312 

1.311 

1.316 

1.317 

1/20 

.36  H 

.62  H 

1.24  H 

2.2  H 

8.64  H 

.340 

.286 

.417 

.283 

.304 

3.647 

3.787 

3.439 

3.795 

3.742 

1.305 

1.326 

1.276 

1,327 

1.319 

Tabu  lb.  Same  as  Table  I  a  except  d  1.5.  The  symbol  S  indicates 
a  separation  transition.  A  hlank  entry  indicates  no  transition  (the 
upstream  tlow  is  not  subcritical). 


R 

Ho 

1 

2 

4 

8 

32 

1/4 

_ 

1.88  H 

3.14  S 

5.6  S 

22.7  H 

.355 

.036 

.016 

.245 

1.358 

1.603 

1.619 

1 .439 

1.167 

1.297 

1.306 

1.209 

1/6 

.91  H 

1.51  H 

2.60  H 

4.80  H 

27.8  H 

.703 

.2435 

.176 

.179 

.847 

1.296 

1.756 

1.823 

1.820 

1.152 

1.105 

1.297 

1.328 

1.327 

1.052 

1/10 

.70  H 

1.16  H 

2.12  H 

4.00  H 

15.0  H 

.51 

.31 

.35 

.344 

.314 

1.875 

2.165 

2.107 

2.119 

2.16 

1.236 

1.332 

1.312 

1.316 

1.33 

1/20 

.48  H 

.82  H 

1.48  H 

2.88  H 

10.8  H 

.344 

.358 

.338 

.37! 

.342 

2.975 

2.945 

2.998 

2.918 

2.980 

1.328 

1.322 

1.333 

1.316 

1.329 

It  will  now  be  shown  that  neither  Lx  nor  L2  can 
change  sign  farther  downstream  in  the  region  where 
their  amplitudes  are  finite.  Suppose  the  contrary  were 
true,  and  let  x’  denote  the  smallest  x  at  which  either 
Li  or  L2  vanish.  If  L2(x')  =  0  with  /.,(.*')  >  0  then 
(3.1)  yields  V  =  1,  and  (3.2)  yields  D  =  W0  +  Lx 
>  IV0,  which  contradicts  the  initial  assumption  (Fig. 
2)  that  D( x)  W'o-  If,  on  the  other  hand,  Li(x')  =  0 
with  L2(x')  >  0  then  (3.3)  gives  U  ~  1,  whereas  (3.1 ) 
gives  U  >  1,  and  thus  a  contradiction  is  again  obtained. 
Therefore  all  admissable  solutions  of  the  nonlinear 
equations  ( 3. 1 )— ( 3.3 >  have  positive  Lu  L2  if  (3.7)  is 
satisfied. 

The  solution  of  the  nonlinear  equations  is  obtained 
by  multiplying  (3.1 )  with  1  -  L2,  by  multiplying  (3.3) 
with  ( 1  -  L2)3U2,  and  by  eliminating  both  LtU  and 
U  from  (3.3)  to  obtain  a  seventh-order  polynomial 
equation  for  the  single  unknown  L2.  The  very  large 
number  of  terms  in  the  coefficients  of  this  polynomial 
were  assembled  by  a  symbolic  computer  program,  and 
then  the  L2  root  was  computed  by  a  careful  numerical 
search  for  the  zero-crossing  point  of  the  polynomial. 
Since  this  starts  from  a  point  far  upstream  where  L2  is 
very  small,  very  small  L2  increments  must  be  taken  in 
the  search,  and  then  small  D  increments  must  be  taken 
to  remain  on  the  same  branch  of  the  seventh-order 
polynomial  as  one  proceeds  downstream.  The  run  ter¬ 
minates  at  that  D  =  Z)cn,  for  which  there  is  either  no 
root  (“hydraulic  transition"  point)  or  one  of  the  ve¬ 
locities  (Ui)  turns  negative  (“separation”).  Tables 
la.b.c  list  the  solution  for  various  values  of  |3,  R ,  WQ, 
and  the  letter  H  or  S  indicates  the  type  of  critical  point. 
For  example,  when  /3  =  1.5  (Table  lb),  R  =  l/io,  W'0 
=  1 ,  a  hydraulic  transition  ( H )  occurs  ( first )  when  the 
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Tabi  f.  lc.  Same  as  Tables  la,  lb  except  J  --  1.25. 


Ho 

R 

1 

2 

4 

8 

32 

1/4 

1/6 

— 

— 

— 

6.08  S 

21.6  S 

.027 

.007 

1.527 

1.543 

1.271 

1.280 

1/10 

.94  H 

1.59  H 

2.68  H 

4.88  H 

18.24  H 

.758 

.336 

.174 

.179 

.198 

1.241 

1.663 

1.825 

1.820 

1.801 

1.084 

1.256 

1.329 

1.327 

1.318 

1/20 

.67  H 

1.12  H 

1.96  H 

3.68  H 

14.0  H 

.546 

.410 

.336 

.337 

.347 

1.981 

2.209 

2.327 

2.326 

2.311 

1.230 

1.295 

1.331 

1.331 

1.326 

minimum  channel  width  is  D  -  0.70,  and  at  this  section 
{/,  =  0.51,  U  =  1.875,  U2  =  1.236.  The  accuracy  of 
the  roots  was  checked  by  using  those  values  to  compute 
the  local  volume  flux,  and  by  comparing  the  result  with 
the  flux  at  x  =  —  oo. 

The  downstream  evolution  of  the  jet  is  illustrated  in 
Fig.  5a  for  WQ  -  8,  R  ~  1/6,  0  =  1.5.  This  shows  that 
as  x  increases  from  -oo  (where  D  =  8),  »ic  value  of 
U\  decreases  and  a  hydraulic  transition  point  occurs 
at  D  =  4.80  with  minl/i  >  0.  At  a  larger  R  (Fig.  5b), 
with  other  things  being  equal,  a  separation  (U,  =0) 
transition  occurs  first  at  D  =  5.6.  In  both  cases  a  jet 
( U  >  U2>  Ut )  forms,  but  Tables  la,b,c  indicate  that 
the  separation  regime  is  relatively  restricted,  occurring 
only  for  the  larger  Rossby  numbers  and  for  the  smaller 
topographic  effect  (0  ~  1 ) .  The  table  also  shows  that 
decreasing  R,  with  fixed  Hq  and  0,  requires  a  decreasing 
Z>crit  ( i  e.,  a  greater  strait  convergence)  for  a  critical 
state,  and  the  accompanying  maximum  jet  velocity  U 
is  increased. 

When  blocking  occurs  in  Fig.  2  the  new  state  that 
arises  could  be  obtained  by  generalizing  the  procedure 
in  section  2,  but  the  algebra  then  becomes  formidable. 

Interesting  variations  in  Fig.  2  and  Fig.  4  might  be 
obtained  by  less  restrictive  shelf  widths  (IF,,  W2),  in 
which  case  separation  might  be  more  prevalent.  The 
effect  of  a  sill  might  be  obtained  by  allowing  0  to  vary 
with  D.  It  may  also  be  possible  to  replace  the  rigid 
bottom  in  the  central  channel  ( Fig.  2 )  by  an  interface 
beneath  which  the  fluid  density  increases,  and  above 
which  the  potential  vorticity  is  (piecewise)  uniform. 
Then  the  important  role  of  buoyancy  on  the  blocking 
could  be  ascertained. 

4.  Conclusions 

A  jet  can  be  generated  by  inviscid  cross-isobath  flow 
through  the  converging  half  of  a  strait,  and  two  types 
of  transition  may  allow  the  generated  vorticity  to  ap¬ 


pear  in  the  diverging  half.  The  local  flow  may  become 
hydraulically  critical  in  which  case  the  generated  to¬ 
pographic  blocking  wave  modifies  the  upstream  state 
as  well  as  the  downstream  state  (Fig.  4  and  section  2 ). 
In  the  second  kind  of  transition  an  everywhere  sub- 
critical  jet  may  separate  because  of  the  relative  vorticity 
generated  in  the  narrow  part  of  the  strait  [cf.  Fig.  5b 
and  Stern  and  Whitehead  ( 1990)]. 

Both  transitions  are  examples  of  purely  inertial  (as 
contrasted  with  diffusive)  mechanisms  for  generating 
inshore  vorticity  in  the  passage  of  currents  through 
straits.  They  may  therefore  be  relevant  to  the  expla¬ 
nation  of  the  large  cyclonic  vorticity  in  the  Straits  of 
Florida,  and  to  the  explanation  of  the  reformation  of 
the  western  boundary  jet  in  the  vicinity  of  the  Straits 
of  Yucatan. 


D 


D 


Fig.  5.  (a)  Downstream  variation  of  velocity  in  Fig.  2  when 
min£)(x)  equals  the  critical  value  for  a  hydraulic  transition,  and 
when  W0  =  8  =  £>(-  x  ),  R  =  i/g,  f}  -  1.5.  I/,  is  the  uniform  velocity 
on  the  upper  shelf.  Ui  is  the  velocity  on  the  lower  shelf,  and  V  is  the 
uniform  velocity  between  the  two  potential  vorticity  interfaces.  A  jet 
with  U  >  t/j,  V  >  V,  forms  at  the  transition  point  D  -  4.8.  (b) 
Same  as  (a)  except  that  R  =  1/4.  In  this  case  a  separation  transition 
(C'i  =  0)  at  D  ~  5.6  occurs  before  a  hydraulic  transition. 


854 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Voi  (  :mk  23 


Similar  kinds  of  transitions  could  occur  in  a  purely 
continental  type  of  topography  (cf.  Stern  1991,  Fig. 
2),  such  as  is  obtained  by  removing  the  sidewall  in 
deep  water  and  curving  the  sidewall  on  the  shelf.  1 
would  like  to  take  this  opportunity  to  correct  an  omis¬ 
sion  in  that  paper  (Stern  1991).  A  term  r  sinhAe 
+  coshAe  should  appear  in  the  denominator  of  the  in¬ 
tegral  of  Eq.  (4.9),  but  this  term  equals  unity  in  the 
long-wave  (A  -*>  0)  limit,  and  consequently  the  omis¬ 
sion  affects  none  of  the  equations  following  (4.9  ). 
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APPENDIX 

A  Smooth  Continental  Slope 

For  previously  mentioned  reasons  it  is  important  to 
compare  the  critical  blocking  condition  for  the  escarp¬ 
ment  model  in  Fig.  3a  with  the  value  obtained  for  a 
smoother  nondimensional  cross-stream  topography 
h{y)  that  is  also  independent  ofx.  For  this  calculation 
it  is  convenient  to  take  the  distance  Ws  =  1  of  the 
upper  wall  from  the  y  =  0  datum  (defined  below)  as 
the  unit  of  length,  with  _v  =  - D(x )  as  the  ordinate  of 
the  curved  lower  wall,  and  with  D(- oo )  =  U'o  as  a  free 
parameter.  The  uniform  velocity  at  x  -  -oc  is  again 
taken  as  the  velocity  unit,  that  is,  U0  -  1 . 

If  8  =  0{y)  denotes  the  ordinate  at  x  =  -oo  of  the 
steady  streamline  whose  ordinate  is  y  at  D(x),  then 
the  Lagrangian  volume  conservation  equation  in  non- 
dimensional  units  is 


,  ,  d$ 

«(>’)  =  «T, 
dy 

(A.1) 

h{8) 

h(y)’ 

(A.2) 

and  the  potential  vorticity  conservation  equation  is 

dy  R  W  )  ’ 

(A. 3) 

where 

R-  Uo 

fw\  ’ 

( A.4) 

The  solution  of  ( A.  1 )-( A. 3 )  for  y  = 
the  boundary  conditions 

y{  8)  must  satisfy 

F(l)=  1, 

(A. 5) 

y(-W0)  =  ~D(x), 

(A.6) 

for  a  given  W0  and  for  each  D  =  D{x).  The  numerical 
solution  is  greatly  facilitated  by  replacing  (A.6)  with 
an  assigned  value  of  u(  1 ).  and  then  “marching”  to  8 
=  -  fV0,  at  which  point  the  computed  y  gives  the  value 
of  D  for  each  w(  I ).  The  calculation  proceeds  with  suc¬ 
cessively  smaller  values  of  u(  1 )  <  1  until  the  minimum 


Ta Bt  I  2.  Minimum  channel  width  U>  />,,,)  for  bh'cking  as  a 
function  of  the  inverse  continental  slope  tv).  Here  lit  D)  is  the  cor¬ 
responding  velocity  at  the  curved  (lower)  wall,  and  id  I  i  is  the  velocity 
at  the  upper  wall:  R  --  I /?<).  3  -  1.5.  H'„  ~  I.  See  Appendix  for 
definitions 


S 

Aw, 

u(  /)) 

X 

.50 

.42 

3.3 

10 

.499 

.400 

3.29 

5 

.513 

222 

3.30 

2.5 

.535 

.286 

2.71 

1 

.687 

.166 

1.43 

D  [as  a  function  of  «( 1 )]  is  reached.  This  gives  the 
critical  amplitude  of  the  obstacle  necessary  to  block 
the  upstream  flow,  and  to  cause  a  hydraulic  transition. 

Calculations  were  made  for  the  hyperbolic  tangent 
profile 

h(y)  =  ~P  tanh(sy)  +  p  tanhi  +  1,  (A. 7 ) 

<3-1 

P  tanh(sH-’o)  +  tanhv  ’ 

whose  maximum  slope  defines  the  y  =  0  datum  level, 
and  which  satisfies  h(  1)  =  1.  h{-WQ)  =  <3.  The  free 
parameter  s  is  a  measure  of  the  inverse  width  of  the 
continental  slope,  and  s  -*•  oo  gives  the  escarpment 
limit  in  Fig.  3a. 

A  second-order  Runge-Kutta  integration  was  used 
to  solve  Eqs.  (A.l  )~(A.3),  and  the  step  size  used  was 
checked  by  comparing  the  total  volume  transport  at  a 
with  an  analytic  calculation  at  x  =  -oc.  Table  2  gives 
the  results  for  four  finite  s  when  W0  -  1 ,  R  =  1/30,  <3 
=  1.5.  The  s  =  00  entry  was  obtained  from  (2.4b)  for 
the  escarpment  model  in  section  2,  and  we  conclude 
that  the  critical  value  of  D  for  blocking  is  essentially 
the  same  in  the  interval  5  <  s  <  oo ,  thereby  validating 
the  asymptotic  consistency  of  the  escarpment  model. 
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